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INTRODUCTION 

Since  the  advent  of  nucl.ear  weapons,  the  investigation  of 
the  ground  shock  effects  from  these  devices  has  been  the  subject  of 
much  research.  Reliable  methods  for  estimating  the  stresses  and  the 
ground  motions  transmitted  from  a  postulated  nuclear  detonation 
through  earth  materials  are  necessary  for  the  cost  effective  design 
of  hardened  strategic  systems. 

Methods  currently  employed,  for  the  evaluation  of  ground 
shock  effects  are  generally  inadequate  from  the  standpoint  of  the 
system  designer.  Computer  code  simulations  arc  designed  to  predict 
the  history  of  ground  motion  from  the  initial  moment  of  the  explo¬ 
sion,  thus  requiring  an  elaborate  description  of  the  characteristics 
of  the  material  under  high  pressure  and  temperature.  Because  of 
their  complexity,  the  extension  of  a  numerical  calculation  into  the 
range  beyond  the  closc-in  field  is  costly  and  often  unreliable  due 
to  the  accumulating  errors  of  the  approximation.  Empirically  deter¬ 
mined  formulae,  derived  from  scaled  high  explosive  and  past  nuclear 
test  events,  are  more  generally  used  to  provide  quantitative  de¬ 
scription  for  specific  weapons  effects  conditions;  however,  these 
formulae  often  cannot  be  extended  to  account  for  varying  geometry  of 
bursts  and  changes  in  the  earth  properties. 

This  situation  is  further  complicated  by  the  fact  that  the 
shock  wave  from  a  surface  burst  preserves  the  close-in  characteris¬ 
tics  of  the  individual  explosion  throughout  the  entire  range  of  in¬ 
terest.  Thus,  slight  deviations  in  the  input  values  from  the  ideal¬ 
ization  of  the  real  environment  can  result  in  seemingly  unrealistic 
results.  The  sensitivity  of  the  solution  to  variations  in  the  source 
description  can  be  mathematically  related  to  a  singularity  of  the 
solution  at  the  origin. 

This  paper  describes  an  attempt  to  improve  the  situation  by 
utilizing  a  simple  model  which  is  adequate  to  describe  the  salient 
features  of , the  ground  motion  history  outside  the  close-in  range.  It 
is  essential,  to  this  purpose,  to  recognize  the  fact  that  the  elastic 
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solution  cannot  account  for  flie  behavior  near  the  explosion  source 
and  that  the  close-in  effects  are  taken  into  account  by  an  appropri¬ 
ate  fit  of  the  singularity  of  the  solution.  This  approach  is  justi¬ 
fied  when  these  effects  are  confined  to  a  small  region  in  space  and 
thus  provides  the  key  to  the  analysis. 


SOLUTION  FOR  THE  CONTAINED  BURST 


Ground  shock  effects  research  has  centered  around  the  fully 
contained  underground  weapons  testing  because  of  limitations  on  test¬ 
ing  in  the  atmosphere.  The  analytical,  solution  for  the  contained  ex¬ 
plosion  is  studied  to  provide  a  simple  check  case  for  verification 
of  the  elastic  model  and  to  determine  the  experimental  scale  factors 
required  for  fitting  the  close-in  effects.  Since  the  particle  motion 
is  most  frequently  measured  in  underground  experiments,  the  particle 
velocity  parameter  will  be  used  in  this  study. 

The  general  solution  for  spherically  diverging  waves  in  an 
elastic  medium  expressed  in  terms  of  the  particle  velocity  u,r,t) 
is  (1) 

u(r,t)  =  -ry  [F(r  -  ct)  -  rF'(r  -  ct )]  (l) 


where  r  is  the  distance  from  the  origin,  t  is  time,  and  c  is  the 
dilatations!  wave  velocity.  F(r  -  ct)  is  an  arbitrary  function  of  the 
argument  (r  -  ct).  A  prime  denotes  a  derivative  with  respect  to  the 
argument  and  a  dot  denotes  a  derivative  with  respect  to  time. 

Applications  of  this  solution  to  waves  generated  by  explo¬ 
sions  were  studied  by  a  variety  of  authors  (1,2).  The  explosion  is 
commonly  modeled  by  applying  a  pressure  to  the  surface  of  a  spherical 
cavity.  The  results  are  generally  found  to  be  inadequate  for  de¬ 
scribing  weapons  effects  because  the  far  field  solution  predicts 
spatial  peak  attenuations  of  nearly  r~l  whereas  measurements  show- 
rates  of  about  r~2  . 

A  close  examination  of  Equation  1  shows  that  by  simply 
specifying  the  arbitrary  function  F(r  -  ct)  in  a  proper  form,  the 
•first  terra  (r~2  terra)  of  the  solution  can  be  made  to  dominate  and 
thus  provide  an  expression  that  at  least  describes  the  peak  values 
of  an  experiment  correctly.  The  function  F(r  -  ct)  was  chosen  as 


F(r  -  ct)  -  A.e 


[2) 


because  of  its  simple  form  and  fit 
expression  for  particle  velocity  is 


u(r,t)  = 
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measurements.  The  final 


(3) 


where  A  and  a  are  quantities  to  be  determined  that  are  charac¬ 
teristic  of  the  source.  The  peak  particle  velocity  occurs  at  r  =  0 
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and  attenuates  as  r~  '  which  is  con si :: Lent  villi  measurements.  The 
wave  form  is  an  exponential  which  arrona  with  the  typical  pulse  de¬ 
termined  by  weapons  testing. 

Conventional  cube  root  scaling  is  introduced  into  Equa¬ 
tion  3  without  (‘hanging  its  dimensions  in  order  to  determine  the 
variation  with  the  weapon  yield  W  .  Equation  3  becomes 
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where 


A  -  A  W 
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V3 


r  W1/3  ,  T  =  t  W1/3  ,  a  -  or  W'  ’  ' 3 
O  O'  o 


The  weapon  parameters  A  and  cy  arc  easily  determined 
from  the  scaled  field  measurements 'by  noting 
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(5) 


where  una>:  is  the  maximum  partic.ic  velocity  and  t4  is  the  posi¬ 
tive  duration  of  a  measured  wave  at  the  sealed  range  rQ  ~  a0  . 
Evaluations  of  A0  and  a*o  show  that  they  are  constants  not  de¬ 
pending  on  test  material  or  weapon  yield. 

A  comparison  between  the  wave  form  calculated  by  Equation  4 
and  a  wave  form  from  the  Salmon  (5  kt)  (3)  event  in  a  salt  medium  is 
shown  in  Figure  1.  The  elastic  calculation  necessarily  has  an  abrupt- 
rise  on  the  initial  portion  because  of  the  need  to  make  F(r  -  ct) 
dominate  its  derivative  in  the  neighborhood  of  the  peak  value.  The 
peak  particle  velocity  calculated  using  constant  A0  is  compared 
with  several  nuclear  bursts  in  a  variety  of  earth  materials  in 
Figure  2. 


SOLUTION  FOR  THE  SURFACE  BURST 

The  object  of  this  section  is  to  utilize  the  classical 
lineal'  clastic  model  for  the  prediction  of  the  characteristics  of 
the  ground  shock  waves  generated  by  explosions  at  the  surface  of  the 
earth.  The  use  of  this  model  for  this  purpose  has  not  been  fully 
explored  because  of  the  lack  of  an  adequate  analytical  solution  to  be 
utilised.  Since  the  mathematical  problem  of  the  explicit  solution  of 
the  surface  burst  explosion  is  very  difficult,  even  for  the  simplest 
model  of  linear  elasticity,  the  existing  analytical  solutions  to  this 
problem  are  either  incomplete  or  unrealistic  for  all  practical  pur¬ 
poses;  thus,  a  more  realistic  and  explicit  solution  had  to  be  devel¬ 
oped.  It  is  essential,  to  this  purpose,  to  recognize  the  fact  that 
the  solution  of  the  surface  burst  problem  is  very  sensitive  to  the 
input  conditions.  Mathematically,  this  sensitivity  is  related  to 
the  singularity  of  the  solution  at  the  explosion  source,  and  the 
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solution,  to  be  useful,  should  be  constructed  such  that  the  singu¬ 
larity  properly  represents  the  characteristics  of  the  source. 

Consider  the  displacement  potentials  cp  and  •!/  determined 
from  the  wave  equations 


■with  propagation  velocities 

.2 

shear  waves  and  A  -  — -  -!- 
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r  or 
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of  the  compressional 
using  the  cylindrical 


(6) 


and 

coor¬ 


dinate  system  r  and  z  with  axial  symmetry.  The  formal  solution 
of  Equation  6  to  the  surface  burst  problem  can  be  obtained  in  the 
form  of  the  double  integral  of  the  Laplace -Hansel  transforms.  The 
direct  derivation  of  the  explicit  expression  of  cp  and  from 
this  formal  solution  is  difficult  for  the  general  region  of  r  >  o  ^ 
but  it  can  be  obtained  on  the  axis  r  -  o  in  the  following  asymp¬ 
totic  expansion  form: 


<P 


> 


(7) 


where  a.  and  b.  (i  =  1,2,...)  are  functions  of  t  -  ~  and  t  -  " 
1  i  ;  c  c 

s 


respectively  and  are  determined  from  the  characteristics  of  an  indi¬ 
vidual  input  source  condition.  Equation  7  exhibits  the  nature  of  the 
singularities  on  the  input  source. 

Approximate  expressions  for  cp  and.  ^  for  r  >  o  arc  . 
constructed  from  the  sum  of  the  elementary  solutions  of  the  wave 
equation  which  have  singularities  at  the  origin  that  arc  matched  to 
the  solution  (Equation  7)  above  are  / 
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where  cpj_  and  \J-j_  (5  =  1,2,...)  are  determined  from  and  b-^ 

(i  =  1,2,...)  by  comparing  the  terms  with  the  same  powers  in  z  at 
r  =  o  with  those  in  Equation  7.  The  horizontal  and  vertical  dis¬ 
placements  u  and  v  are  then  determined  by 
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Tne  accuracy  of  the  approximation  depends  upon  how  many 
terms  are  retained  in  Equation  8.  The  result  of  the  application  of 
this  method  to  the  simpler  case  of  the  liquid  half-space  with  only 
one  wave  shows  good  agreements  with  experimental  data  is  attained 
with  two  terms  of  the  approximation  formula  (ii).  a  cursory  examina¬ 
tion  of  the  elastic  case  with  two  terms  of  the  approximation  is  bein^ 
made  and  compared  with  test  data  for  rock  environments.  Tile  results'3 
of  this  correlation  are  shown  in  Figure  3,  where  the  horizontal  ac- 
cel orations  u  near  the  surface  are  compared  with  the  MINE  ORE  (5N 
high  explosive  test  data. 

CONCLUSIONS 

General  formulae  for  the  ground  motions  from  an  explosion 
source  were  derived  and  their  results  compared  with,  high  explosive 
and  nuclear  weapons  test  data.  The  compari  sons  show  reasonably  good 
agreement  may  be  obtained  by  using  a  simple  elastic  description  for 
the  earth  as  long  as  the  explosion  is  properly  modeled. 
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Fig.  1.  Comparison  of  calculated  particle  velocity-time  history 
with  test  data  from  Event  SAJl-iON  at  r  -  166  meters 
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Fig.  2.  Comparison  of  calculated  peak  particle  velocity  with 
test  data  from  nuclear  events  in  rock 
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Fig.  3-  Compari son  of  calculated  horizontal  acceleration-time 
history  with  test  data  from  MIKE  ORE  Event,  MINE  SHAFT  Scries 


9 


